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Abstract. We show that a correct formulation of the cold collision frequency shift for two 
photon spectroscopy of Bose-condensed and cold non-Bose-condensed hydrogen is consistent 
with experimental data. Our treatment includes transport and inhomogeneity into the theory 
of a non-condensed gas, which causes substantial changes in the cold collision frequency shift 
for the ordinary thermal gas, as a result of the very high frequency (3.9kHz) of transverse 
trap mode. For the condensed gas, we find substantial corrections arise from the inclusion 
of quasiparticles, whose number is very large because of the very low frequency (10.2Hz) of 
the longitudinal trap mode. These two effects together account for the apparent absence of a 
"factor of two" between the two possibilities. 

Our treatment considers only the Doppler-free measurements, but could be extended to 
Doppler-sensitive measurements. For Bose-condensed hydrogen, we predict a characteristic 
"foot" extending into higher detunings than can arise from the condensate alone, as a result of 
a correct treatment of the statistics of thermal quasiparticles. 



1. Introduction 

Bose-Einstein condensation of spin-polarized hydrogen, achieved in 1 998 |l^ ^ ||, ||| as the 
culmination of a quarter of a century of development of cooling, trapping and spectroscopic 
techniques, still presents certain puzzling features. The most notable of these occurs in the 
spectroscopic method used to detect the condensate, based on two-photon excitation from the 
IS to the 2S state of the hydrogen atom in the presence of a significant background density 
of hydrogen atoms — a density which of course becomes very large when the hydrogen Bose 
condenses. The cold collision frequency shift — a shift of the transition frequency proportional 
to the density of the gas — is measurable even when the gas is merely dense, but not Bose- 
condensed, and is very significant for the condensate. This has led to the use — suggested by 
Shlyapnikov (see reference [19] of [|l]]) — of the cold collision frequency shift as a method 
of measuring the proportion of atoms which occur at a given density in the system of cold 
hydrogen atoms. The magnitude of the shift is expected to be twice as much for a given 
density for the non-condensed gas as for the Bose-condensate, as a result of the different 
coherence properties of the two systems. However, the frequency shifts observed have led the 
experimenters to the conclusion : 

... if we had assumed that the frequency shift in the condensate is only half as large 
as for a normal gas at the same density, as one would expect for a condensate in 
a single quantum state, the density extracted from the spectrum would be twice 
as large. This would imply that Nc ~ 6 x 10^ and yield an unreasonably high 
condensate fraction of 25%. 
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The apparent absence of the "factor of two" has led to the speculation that the hydrogen 
condensate is not in fact a true condensate, but is in some unexplained way incoherent, so 
that the behaviour is qualitatively similar to that of a non-condensed system, and Cote and 
Kharchenko [|^] have given a more exotic explanation, based on droplet formation at the edge 
of the condensate. Our aim in this paper is to take the view that the hydrogen Bose-Einstein 
condensate has no exotic features, but to make a careful investigation of the physics of the 
two-photon excitation process used to probe the system, both when it is condensed and when 
it is non-condensed. In our view, the apparent absence of the "factor of two" implies that 
something is wrong with the theoretical treatment of the cold collision frequency shift in 
either the condensed system or the non-condensed system or possibly in both, as they have 
been so far formulated. 

We have therefore reviewed the conventional theory in Sect. 1, and have shown that for 
a non-condensed system it is assumed that the system is homogeneous. The treatment of the 
condensed system [H, || does not make this assumption, but requires the system to be a pure 
condensate. Neither of these requirements can be considered to be adequately met in current 
experiments, since the system is tightly trapped and relatively warm. 

In Sect.2 we develop a theory of the cold collision frequency shift in inhomogeneous 
non-condensed systems, and show that in the published experiments [ [l0| , ^] the effect of 
the inhomogeneity will be very significant, and will tend to reduce the "factor of two" 
significantly. The main reason is quite straightforward. During the time of the exciting pulse 
the atoms move through a wide range of densities, and the coherence between the excited 
state and the ground state, necessary to produce the "factor of two", cannot be established. 

In Sect. 3 we turn to the condensed system, and develop a formalism which uses a 
quasiparticle basis to give some linear partial differential equations for the excitation as a 
function of space and time. We use the very anisotropic aspect ratio of the hydrogen trap — 
about 400:1 — to simplify our formalism to a relatively simple set of formulas from which 
quantitative predictions can be made. Our results confirm those of Killian for a pure 
condensate, but also give results which are valid when a significant quasiparticle component 
is present. In essence we show that the quasiparticle component can be up to 20% of the 
apparent condensate occupations, and that the condensate and the quasiparticles give rise to 
independent and qualitatively different signals. For a given excitation frequency, the two 
signals are produced from different locations in the system — this can be seen as arising 
because there is not only a cold collision frequency shift but as well a splitting which gives 
rise to two frequencies, the two components corresponding to the two signals. The situation 
is illustrated in Fig.|[ In fact, at the highest detunings, the signal arises only from the 
quasiparticle component, and gives rise to a characteristic "foot" protruding into the high 
frequencies above an otherwise smooth spectrum, and indeed there is some experimental 
indication of such a feature — see Fig.^. 

We conclude that if the our calculations for condensed and non-condensed experiments 
are taken into account that the data of appear to be consistent with the following: 

i) The calculated value of the 1S-2S cross section. 

ii) The condensate is coherent, and is accompanied by the appropriate distribution of 
quasiparticles for the measured temperature. 

iii) The non-condensed system is incoherent, with 52 = 2. 
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2. The cold collision frequency shift for a homogeneous system 

In an inhomogeneous and possibly Bose condensed system, we consider the excitation by a 
coherent process. One photon excitation by a coherent microwave field is usual for an atomic 
clock, but in the case of a hydrogen Bose-Einstein condensate two-photon excitation is used to 
excite from the IS to the 2S level. The exciting field is pulsed for a time very long compared 
to the transition frequency, but so short that only an infinitesimal proportion of the atoms are 
excited to the upper level. 

2.1. Second-quantized Hamiltonian 

For the case under consideration the second-quantized Hamiltonian can be written in the form 

j d^yi (T + y(x) + fiLJi) V^i + iPI (T + F(x) + huj2) ^2 



u 



(1) 

i) The labels 1, 2 represent the IS and 2S states of hydrogen, and huji, fiuj2 are IS and 2S 
energy levels used in the experiment. 

ii) The kinetic energy operator is T = — V^/2m 

iii) This form assumes the trapping potentials have the same form ^(x) for the IS and the 
2S states of hydrogen. 

iv) Since the process being considered is a two-photon process, the driving term is 
proportional to the square of the laser field, which is assumed to have frequency 57/47r, 
and wavenumber k/2, so that we can write the driving term using 

£;(x) = £;(! + cos(k- x)). (2) 

The spatially independent term gives the Doppler-free excitation, while the cosine term 
gives the Doppler-sensitive excitation. In the experiments both forms of excitation are 
used — however in this paper we shall treat mainly the Doppler-free case, for which we 
have been able to develop a relatively simple formalism. 

v) The interaction coefficients are given in terms of scattering lengths by 

u = Anais-ish^ /ni (3) 
V ^ 47rais-2S?i^/w (4) 
w = 47ra2s-2S?i^/w (5) 

The scattering lengths have been calculated in ^ to be 

ais-is = 0.0648nm (6) 
ais-2S = - 2.3mn, (7) 

and the experimental value (calculated assuming that the theory of the cold collision 

frequency shift in a non-condensed trapped system is the same as that of a homogeneous 

trapped system) 

ais-2S = - 1.4±0.3nm. (8) 
is reported in [ [lo| ] . 

vi) The coefficient w will be set equal to zero, since its effect is in practice negligible with 
the densities of 2S atoms appropriate for the experiment. 
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2.2. Background 

The approach used by Kilhan || considers an initial many-body wavefunction composed 
as a symmetrized product of one-body wavefunctions \ IS). The resuh of the laser excitation 
is to rotate the spin wavefunction slightly into the 2S subspace, so the one-body wavefunction 
is transformed 

\1S) ^t\lS)+r\2S), (9) 

with |tp + |rp = 1, and r <C 1. In field theoretic language (in the Heisenberg picture) the 

initial state is described by field operators ■01 (x), V'2(x) with (-i/'2(x)'02(x)) ~ 0, and these 
transform to 

V'i(x)-> tM^)+rM^) (10) 

V'2(x)-> -rM^)+tM^)- (11) 

and the population densities ni(x) = (-i/'|(x)-0i(x)) transform as 

ni(x)^ni(x)~57i(x) =|ipni(x), (12) 

n2(x)^5n(x) =|rpni(x). (13) 

The mean energy after this process is given by the sum of kinetic, potential and self energy 
terms. Since the trapping potentials for the two states are assumed the same, there is no change 
of trapping energy, and if the excitation is not strongly dependent on space there is no change 
in the kinetic energy term. Under these conditions the main energy change is given by the 
interaction energy, which changes by an average amount 

(1^1' - 1) + v\tW) (V'I(x)0l(x)V'i(x)V'i(x)) « 52(x,x)(« - u)n,{^)Sn{^). (14) 

Here we have used the spatial density correlation function 

ni(x)ni(x') ■ ^ 

If the initial quantum state is noncondensed, it can be written in terms of a set of orthonormal 
one particle wavefunctions 0r-(x) with occupations Nr which are either or 1, so that 

ni(x) =Y,\M^)\^Nr, (16) 

r 

g2(x,x)ni(x)2 = 2j2\M^)\^\M^)\^NrNs (17) 

«2m(x)2, (18) 

where the approximate final result is valid provided the occupation is spread over very 
many modes, so that the missing term with r = s is negligible. Thus we find in this case 
32 (x, x) 2. This result is also obviously true for an ensemble of noncondensed quantum 
states, such as a thermal noncondensed state. 

If the state is condensed in the sense that only the quantum state r = rg is occupied, and 
^ 1 , then the same calculation leads to .92 (x, x) 1, so that the energy shift is only 
half of the value expected from a noncondensed state — the "factor of two" thus does appear 
when the trapping and kinetic parts of the Hamiltonian can be neglected. 
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2.2.7. Issues for a realistic condensate There are three main issues which need to be 
addressed 

i) The energy shift is proportional to ni(x), and is sufficiently large in practice for it to 
be used as a spectroscopic method of measuring the density profile of a very cold and 
possibly condensed cloud of hydrogen atoms. Thus the amplitude t for transition to the 
2S state is of necessity spatially dependent. This issue has been addressed in ^] for the 
case of a pure condensate, but it has been assumed that the results of the homogeneous 
noncondensed system are applicable to the inhomogeneous non-condensed case which 
occurs in the experiments. A full consideration of the effects of inhomogeneity for a 
noncondensed system or a partially condensed system has not been made. 

ii) The temperature of the condensed system is relatively high, so that occupation of the 
lowest quasiparticle levels is significant. We showed previously [||] that this could be 
up to 20% of the apparent condensate population, and that the effective 92 (x, x) in such 
a situation is spatially dependent, varying from a little more than 1 to somewhat larger 
than 2, depending on the position in the condensate. Under these conditions the simple 
argument given above cannot be justified. 

iii) The estimate of the energy shift ( p^ uses an average energy. However, at a given position 
in space, this could be the result of an average of different eigenvalues, and it would be 
these eigenvalues which were selected spectroscopically, not their average, leading to a 
given resonant frequency corresponding to different densities of the gas. 

3. Cold collision frequency shift for an inhomogeneous non-condensed system 

As noted above, the "factor of 2" has been derived only for homogeneous systems, and in the 
experiments ^ homogeneity is not satisfied. One need only note: 

The radial trap frequency is « 3.9kHz 
The axial trap frequency is w lOHz 

The laser pulse length is « 500^s « 2 radial trap periods. 

Thus, during the period of the excitation it will be possible for a significant number of atoms 
to experience a wide range of values of the vapour density, and a treatment must be devised 
which includes the full dynamics of transport in the tightly trapped cloud. We will show that 
a combined transport-excitation equation in 3 position plus 3 momentum variables can be 
derived. 

3.1. Equations of motion 

From the Hamiltonian (|l]) the Heisenberg equations of motion for the field operators are (in 
the Doppler-free case) 

ihiPi [hwi + H{x.)]iji + uVlV'iV'i + «V'2V'iV'2 + (19) 

ini;2 = [riuJ2 + H{x)]ij2 + vijlijiij2 + gEe-'^*^i, (20) 

in which 

H{^)^ _^ + y(x). (21) 

We suppose that the occupation of the IS level is always very much greater than that of the 
2S level, and all of this occupation induced by the application of the coherent driving field 
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E. This means that we can neglect entirely the last two terms in (p^, so that ipi (x, t) can be 
regarded as a known time-dependent operator 
We now define Wigner amplitudes 

/(p,x) = dV(V4(x-y/2)V'2(x + y/2))e-'P-y/^ (22) 

n(p,x) = rfV(V'I(x-y/2)Vi(x + y/2))e-'P-'^/'' (23) 

We wiU also use the momentum integrated quantities 

iV(x)== y"d3pn(p,x)= (^I(x)V'i(x)) (24) 

^^(x) = J d3p/(p,x)= (V't(x)^2(x)) (25) 

which will characterize the system sufficiently for our purposes. We can now derive equations 
of motion using 

i) Hartree-Fock factorization in the form (for example) 

= {4iz)M^')){4{z')M^')) + {4iz)M^')){4{z')M^')), (26) 

with z = X — y/2 and z' = x + y/2. 

ii) Since the terms (-0|(z)-i/'i(z')), (-i/'|(z)?/'2(z')) can only be nonzero for very small 
y = z — z', we also make the approximations of the form 

{ip\{z')Mz')) = n{z') « n(x) + y • Vn(x)/2. (27) 

iii) We can then use the standard Wigner function methods to get the equation of motion 

i^^^^§f^ - H^2 - c^i)/(p, x) + gEer^'^'nip, x) 

+ ih {-(p/m) • + vy(x) • Vp} /(p, x) 
- (2ii - ■!;)A^(x)/(p, x) + vF{x.)n{p, x) 

+ ih (ii + I) ViV(x) • Vp/(p, x) + ih^VF{x) ■ Vpn(p, x). (28) 

3.1.1. The measured signal The measured signal is the total number of atoms excited to the 
2S level, that is 

s{n) = J d^^{4{^)M^))- (29) 

We can use an ansatz 

M^)^ J dVKx,y)^i(y) (30) 

to express the linearity of the process, in the sense of ( p^pT| ), and then we can write 

/(y,x) = {4{y)M^)) ~ [ d'y'r{^,y'){i'l{y)My'))- (31) 



It is reasonable to approximate the correlation function on the right hand side of this equation 
by a locally thermal form 

(V'I(y)V'i(y')) « N exp (- ^^~/^' ) , (32) 
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in which 

L = h/V27rmkT. (33) 

The length L is sufficiently short in a thermal situation to treat the correlation function as 
being essentially local in ( pOpl] ), and thus write 

/(y, x) « (V2^L)3iV(y)r(x, y), (34) 

and using this we can then similarly derive 

s{Q) « / d^x / l/(y.x)p 

^ ' J J (^/2^i)3iV(y) 
5.2. The homogeneous case 

Suppose now that V{x.) = 0, and r7,(p,x) = nhom(p) is independent of x. There is then a 
solution of ( ^8| ) with /(p, x) cx n(p), also independent of x and the corresponding equation 
of motion can be integrated over p to give 

ih— = {n(oj2 - t^i) - 2(w - v)N}F + gENe-'^* (36) 
ot 

Thus in this case we find that 



/(y,x)^fcxp(- ^y^^^^ ) (37) 



and 



where V is the volume of the homogeneous system being observed. We can easily solve the 
equations ( p6| ) to get 

i Vtot(gg)^ sin^{(»- A)r /2} 

272^2 {(f] - A) /2} 
A^tot(gi?)' 



-2TrTS{n-A). (40) 



2(u - v)N 

where A = uj2 — ^^i z (41) 

h 

This confirms the result of ( p4| ) with g2 — 2. 

3.3. One dimensional solutions of the equations of motion 

The equations (^s]) form a six-dimensional set, which are quite formidable to solve 
numerically, so to get an idea of their effect, we have done one-dimensional solutions, with 
trap frequency w,.. The results are presented in Fig|l| The solutions of the equations obtained 
by omitting the "streaming term" on the second line of ( ^8| ) — corresponding to the result 
( |l4l ) — are given in the dashed curves. It can be seen that for the trap frequency frequency of 
= 27r X 3.9kHz , as used in the experiments, the excitation curve is very strongly affected, 
but for the frequency lo'^ ~ uj,./^ ~ 2tt y. 780Hz, one fifth of that frequency, although there 
is a small quantitative difference in excitation, the shape is essentially unchanged. These 
should be compared with the excitation time T = for which LOrT 12^ — 1.95, while 

w'/r/2'K = 0.39. The conclusions we can draw from these one dimensional simulations are 
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Detuning (kHz) Detuning (kHz) 



Figure 1. Excitation as a function of laser detuning. In all cases the laser is turned on for a 
time T = 500/ts. The solid curves are solutions of the equations (ESI), the dashed curves omit 
the "streaming terms" on the second line of (Eg), a) uir = 3.9kHz ; b) u)r = 780Hz. 



limited to the observation that the peak of the soHd curve in Fig.|i]is shifted from the position 
which would be expected if one used the homogeneous system formulation of Sect. 2.2. The 
shift of the peak is such that one would see an apparent "factor of about 1.57", rather than 
the "factor of 2", and this would indicate that the true 1S-2S scattering length should be 
larger than the experiments have actually measured. We are reluctant to treat the value 1.57 
as anything other than a qualitative result, but it seems that this would suggest that it might be 
wiser to use the theoretical value in analysing the condensate data rather than the apparent 
measured value 



4. Cold collision frequency shift for an inhomogeneous condensed or partially 
condensed system 

The treatment in the presence of a condensate takes a completely different form from the 
full thermalized case. As before, the equations cannot be solved exactly, but approximate 
equations of motion which incorporate the features of the quasiparticle spectrum, as 
investigated in [||, can be developed as follows. 



4.1. Approximate equation of motion 



The equations of motion and assumptions used in Sect. 3.1 are valid in this case as well, but 
the Wigner function methodology is not appropriate, since its use depends on hartree-Fock 
factorization. 



4.1.1. Quasiparticle expansions As before, we suppose that the occupation of the IS level is 
always very much greater than that of the 2S level, so that (x, t) can be regarded as a known 
time-dependent operator. However, in this case we represent the time-dependent solution for 
ij^i (x, t) by a quasiparticle expansion of the form 



■!/'i(x, i) = ao(<) 



( iM + ^Y. {«"(0^n(x) + at (i)<(x)} 1 , (42) 

where the time-dependences of the destruction operators are 

ao{t) ^ -i{uji+ ^i/K)ao{t), (43) 

an{t) = -i{en/n)a^{t). (44) 
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The most accurate forms for the ampUtudes u„ (x) , u* (x) , and values of the chemical potential 
yU and quasiparticle energies e„ would be those determined by the gapless formahsm of 
Morgan and coworkers [ pi] , |l2||. However, for simplicity, we will use only the Hartree-Fock- 
Bogoliubov (HFB) method in the Popov approximation, which, as Morgan has shown, is a 
much better approximation than the pure HFB method. Our methodology depends only on the 
existence of approximations of the form (p^[-p4|), and so can be adapted to any approximation 
scheme of this general kind. 

The equation for V'2(x, t) is now linear and inhomogeneous, with operator coefficients. 
We can develop approximate solutions using an ansatz of the form 



ao{t) 



+ terms involving -02 (x, 0), (45) 

where the amplitudes s(x, t), p„(x, t) and qn{^, t) are c-number amplitudes whose equations 
of motion are to be determined. 

This form expresses in a more general sense the concept summarized in (|lo|jllj), that 
the driving field coherently rotates a small proportion of the the field ipi (x, t) linearly into 
the 2S subspace. Here we keep the linearity in the sense that ( p^ is a linear combination of 
components ao{t), ao(^)a„(t), al^aQ{t) in terms of which the field '0i(x, t) can be expressed. 
However, it can be seen that the evolution according to (|o|) depends nonlinearly on the 
IS field, and this leads in the approximation we are using to a dependence on nonoperator 
nonhnear functions of the field ipi, in much the same way as the expression ( |4^ gives 
linearized excitations of the condensate, whose equations of motion nevertheless depend on 
certain averages of nonlinear functions of the field operator. 

4.1.2. Equations of motion From ( p3| ) we can express the coefficients in terms of thermal 
averages 

(aS(t)02(x,O) =^s(x,t), (46) 
(4(i)aUi)V'2(x,<)) - V^iV„p„(x,t), (47) 
(aJ(t)a„(t)V'2(x,t)) = ^{Nn + l)qn{^,t). (48) 
Here we assume the condensate-vapour system in the IS state is in a thermal state, so that 



Using these, it is straightforward to derive the set of coupled equations of motion: 

ihs = {n(L02 -LUi)- ^l + H + v[\^\^ + n]}s + gECe-'^^\ 

+ [Nm[u*jn + V^]PmiNm + l)[Vm + Um]qm^ (50) 
m 

ihpn ^ {h{uJ2 - LUi) - ^ - €n+ H + + n]] Pn + gEUnC"''^^ , 
+ {v[C + S,Vn]} S 

+ W ^ |^m[MnU*„ + Wm^nlPm + (^m + l)[u„W„j + U„W„]g*, | (51) 
m 

ihql = {H^2 -^i)- ^i + en + H + + n]} ql + gi^^e"*"* 

+ {^'K< + r<]}s 

+ W ^ {^™K"m + '"*niK]Pni + {Nrn + l)b>m + UmM*]g*,(x, (52) 
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Here n(x) is the mean noncondensate density 

n(x) = J2 {Nn\un{^)\^ + {N„ + l)|«„(x)|2} . (53) 

71 

In terms of these amplitudes, the occupation of level 2 is 

(4(x,t)Vj2(x,0) = |s(x,t)|2 +^ {7V„|p„(x)|2 + (7V„ + l)|g„(x,0|2} . (54) 

n 

5. Simplified solutions 



The equations (p0[-p2[) are somewhat opaque, so in this section we will develop some 
simplifications so as to see the general structure of their predictions. The simplifications 
we shall make are quite drastic, and thus the results of this section may have only a qualitative 
significance. However, the full equations (^-52) should be rather accurate, and are by no 



means intractable, so that we will always be able to check the approximations explicitly. 
5.1. Development of equations 

In a previous paper [||] we developed approximate wavefunctions for the low lying 
quasiparticle modes, which we expect to dominate that part of the noncondensate fraction 
which exists in the same location as the condensate. It was shown in this paper that the first 
10 eigenfunctions give nearly all the contribution to this part of the noncondensate fraction. 
In this case for such low lying eigenfunctions, we can make some simplifications (which one 
would expect to be generally valid for any kind of quasiparticle formalism): 

i) The occupations Nn are very large, and we can neglect the difference between iV„ and 

Nn + 1. 

ii) The eigenfunctions can be chosen real. 

iii) There is little difference between the two types of eigenfunctions, in the sense that 

w„(x) « -it„(x). 

iv) The quasiparticle energies e„ for these eigenfunctions are negligible compared to the 
other mean field effects. 

Putting in these simplifications, the equations for p„ and —qn become the same as each other, 
so we can write g„(t) ~ — p„(t), and the equations (|5C|-|5^) simplify to 

ihs = {hiuj2 -LOl) - + i?(x) + l'[|C(x)|^ + ?l(x)]} S 

+ .gi?(x)e(x)e-'"* (55) 
ihpn = |ri(w2 - uji) - fi + H{x) + u[|C(x)P + n(x)]|p„ 

+ 4?;^ iV,„u„(x)u„(x)p„, + g£;(x)u„(x)e"'"*. (56) 



5.7.7. Condensate contribution The first of these equations, representing the condensate, 
is now independent of the others. The method of solution is very similar to that used by 
Killian [^, as follows. We replace fi using the Popov and Thomas-Fermi approximations in 
a Hartree-Fock-Bogoliubov theory, from which it follows that the condensate wavefunction 
approximately satisfies 

u|C(x)p = M - "^^(x) - 2un(x) where fi > V{x) + 2un(x); (57) 

— otherwise. 
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Figure 2. Form of the effective potential. 



We can then eliminate /i in the region where the condensate wavefunction is nonzero, which 
is where our formalism is valid and is also where the signal of interest comes from. (Outside 
that region we can make no elimination, so there is an effective potential as in Fig.|— the very 
tight attraction at the centre arises because v is negative, and about 20 times larger than u.) 
Hence we expand in terms of the orthonormal eigenfunctions g{uj,i,x.) of the operator 

Hc^ - ^r— + {v- u)|C(x)|2 + {v~ 2u)n{^). (58) 
2m 

Here Tiuj is the eigenvalue of He, and i represents all the other eigenvalues necessary to 
specify the state. Thus, we write 

s{^,t) = ^s^,i{t)g{u},i,x) (59) 
g£;(x)e(x) = ^/^,,5(c^,i,x). (60) 

iijA 

The solution for s^^i{t) is then clearly 

S^M = /.,.e-^^" — . (61) 

\ rL(UJ2 — UJl + UJ — il) J 

The signal from the experiment is proportional to the total occupation, i.e.. 



2TTt 



_Y^I7- |2 



2 ^|^,»r(5(w2-wi+w-f^). (63) 



LU'l 



The quantity j is given by 

I^a^g j d3x5*(c^,i,x)e(x)i;(x). (64) 

At this stage it becomes important to distinguish between the Doppler free case, where 
£'(x) = f , a constant, and the Doppler sensitive case where £'(x) has a very rapid oscillation 
at the wavelength of the light being used, which makes a simple approximation scheme 
difficult. We do not treat the Doppler-sensitive case in this paper 

5.2. Doppler free excitation 
In this case we write 

/c.,^ y"rf'x.g*(^,i,x)e(x). (65) 
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Figure 3. Plot of the Airy function and its integral 



and follow the reasoning of Killian [|^], that for relatively high energy excitations such as 
occur here, the integral of the product of a relatively smooth function, such as ^ (x) and the 
eigenfunction g*{uj, i, x) is dominated by the contribution at the boundary of the classically 
allowed region, that is from x such that 

nLo = {v- m)|^(x)P + [v - 2u)n{x). (66) 

Qualitatively, this occurs because the eigenfunction g{LLi,i,x.) has a behaviour like that of an 
Airy function near the place where the energy eigenvalue is equal to the potential term in 
Tic, that is, where uj satisfies (^6|). The Airy function and its integral are plotted in Fig.^. 
Combining (^3]) and (^6|), we find that the main contribution comes from positions where 

{v - it)|C(x)|2 + (w - 2u)n(x) 



n 



(67) 



5.2.7. Case of a pure condensate In this case we set n(x) 
as those of Killian ^ ^] for a pure condensate. 



0, and our results are the same 



5.2.2. Quasiparticle contribution The remaining equations can be further simplified by 
noting that because \v\ ^ the spatial variation of p,i(x) will be very much faster than 
that of M,i(x). Using this, we can make the ansatz 

p„(x) =M„(x)P(x) (68) 
and because of the slow variation of Uri(x) we can approximate 

VV(x) « ^i„(x)V2p(x). (69) 
We use this ansatz, and the fact that when 1 is neglected in comparison to Nm 

n{x) fa 2^iV™u™(x)w„(x), (70) 



to conclude that, for any n, the equations ( p6| ) become 
in- 



dt 



{h{uj2 - t^i) - + H{yi) + f [|e(x)|2 + 3n(x)]}p„ + .g£u„(x)e 



(71) 



The fact that we get the same equation for all n means that there are solutions of this equation 
of the form (68) — that is, our requirement that P(x) be independent of n is verified. As in 
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Figure 4. a) Relative condensate (solid line) and quasipailicle (dashed line) contributions to 
the density; b) Detunings generated by condensate (VLq, solid line) and quasiparticles {Q,qp, 
dashed line); see text SectM 



the case for equation for the condensate contribution, we can eliminate the chemical potential 
using define the operator for the quasiparticles 

'Hq= - ^— + {v- 7/)|e(x)|2 + (3i> - 2u)n{^), (72) 
and expand in terms of its eigenfunctions /(tj, i, x), using 

P«(x,i) = ^p„^^^j(t)/(w,i,x) (73) 

g£u„(x) = ^ J„^^^i/([j,i,x). (74) 

We proceed as for the condensate, and find that the quasiparticle signal is 

iin^[(cL;2 — ijJi — 00 — iY)t/ 

— oji + uj ~ ri), (76) 
where 

Jn,uj,i ^ g£ (i^x/*(w,i,x)u„(x). (77) 



2TTt 



In this case, the functions m„(x) are also very smooth compared to the eigenfunctions, and 
the integrals will be dominated by contributions from x such that 

O (.\-., (^;~^)|g(x)p + (3^;-2H)n(x) 

"QP(.xj — tJ2 - t^l H r • (/o) 



6. Experimental Comparison 

The equations ( |66| , [78[ ) can be viewed as the fundamental equations for determining the 
position in the system from which the measured signal (i.e., the total occupation of the 2S 
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Figure 5. Illustration of the contours of excitation expected for a) The quasiparticle component 
and b) The condensate component. The letters A-D correspond to the frequencies f^A-f^D- 
Note that the highest frequency gives no condensate signal. 



State) originates, although they are definitely only approximations to the equations (|50|-p2[). 
It is apparent the the resonance condition for the condensate component, ( |66| ) is significantly 
different from that of the quasiparticle component, (|7^, especially as, v ~ —22u. 

In order to evaluate the full signal from the two components, it will be necessary to 
evaluate the sums ( ^3p^ ) in the case that the eigenfunctions are almost Airy functions, and 
this is not a trivial task. In the case that we can approximate the effective potentials by a 
harmonic approximation, we can show, as in |^ |[], that we can take the intuitive form 



S{Q) cx 



(79) 



in which flc and ilgp are given by ( p7[ ) and (78) 

The measured signal is seen form (^9|) indeed to be the result of an inhomogeneous 
summation — at a given the condensate signal (first term) and the quasiparticle signal 
(second term) come from different regions of the cloud. Using the eigenfunctions from |^] 
we show in Fig. ^ how this happens. Along the line A-A, corresponding to = il^ the two 
signals come from very different positions, while at = Ob there is no condensate signal, 
and the quasiparticle signal is in fact itself derived from multiple positions. What can be seen 
is that the signal with highest detuning is in fact derived from the quasiparticle component of 
the signal. 

The result of this dual source of the excitation signal is to produce a spectrum which 
has two very distinctive components, as shown in Fig.|[ where the condensate component 
is qualitatively the same as one would expect from a Thomas-Fermi treatment, and the 
quasiparticle component is dominated by a large detuning component, which forms a distinct 
"foot", protruding from the from the large detuning end of the spectrum. At a given 
temperature, this feature is more distinct at lower condensate numbers. 

The agreement with the experimental data of [Q] is rather good. We have taken the 
total output signal in the form (79), and used the calculated value, (R), of the 1S-2S 



scattering length, since the results of Sect. 3.3 indicate that the apparent value resulting from 



a measurement in a noncondensed system using the simple g2 = 2 method appropriate to a 
homogeneous system is likely to be significantly lower than the correct value. 



7. Conclusion 



We believe that in this paper we have resolved the apparent anomaly, that the hydrogen 
condensate appeared to be incoherent, and therefore have 172 = 2, rather than 52 = 1, as 
expected for a true condensate. There are two reasons for the appearance of this anomaly: 

i) The use of a formalism only strictly valid for the case of a homogeneous untrapped 
condensate is not justified with the trap geometry employed for the experiment. The 
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Figure 6. Eredifted and measured excitation spectra for a) Parameter regime of the 
experiments M ho||; b) One tentli of tlie number of condensate atoms. Notice tliat tlie "foot" is 
much more prominent in the latter. 



time during which the system is illuminated by the exciting illumination must always be 
significantly shorter than the shortest trap period for this to be valid. 

ii) The more careful treatment of quasiparticle effects leads to higher frequency shifts than 
would be expected from a treatment based on the assumption that the condensate is pure. 

We have also established that the presence of quasiparticles causes a splitting in the resonance, 
rather than a simple shift, and this leads to a a very characteristic "foot" in the excitation 
spectrum, of which there is already some experimental evidence. 

Finally, it must be emphasized that our treatment still makes many approximations. In 
future work we will include mean-field effects in out computations of eigenfunctions, and 
shall also develop methods of solving the full equations so derived, including the Doppler 
sensitive case, which has not been amenable to a simple approximate treatment. 
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